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Abstract 

The quantity Rk = T{K -> eiy)/r{K ^v) studied by the experiment NA62 at CERN 
is known to probe lepton-flavour violating (LFV) parameters of the Minimal Supersymmetric 
Standard Model (MSSM). A non-zero parameter 5^^^ can open the decay channel K — > evr 
and enhance Rk over its Standard-Model value. In the region of the parameter space probed 
by NA62 the contribution from a bino-stau loop diagram is numerically dominant and the mixing 
between left-handed and right-handed staus is important. For large values of the stau mixing angle 
Or the commonly adopted mass insertion approximation is not accurate. We therefore express the 
supersymmetric contribution to Rk in terms of the mass of the lightest stau eigenstate, the mixing 
angle 9^, and other relevant MSSM parameter such as tan/3 and the charged-Higgs boson mass 
Mh and plot the parameter regions constrained by Rk- We further study to which extent Rk 
can be depleted through MSSM contributions interfering destructively with the SM amplitude for 
K — )■ eVe- This lepton-flavour conserving (LFC) mechanism involves the parameter combination 
I'^ii'^jfj?!' which can be constrained with a naturalness consideration for the electron mass or 
with the measurement of the anomalous magnetic moment of the electron. The LFC effect on Rk 
is marginal, an NA62 measurement of Rk significantly below the Standard Model expectation 
would indicate physics beyond the MSSM. 



1 Introduction 

The discovery of neutrino oscillations has shown that individual lepton numbers are not conserved. 
This phenomenon constitutes physics beyond the Standard Model in its original formulation, which 
involves only renormalisable interactions and contains no right-handed neutrino fields. Nevertheless, 
the Standard Model can accommodate neutrino oscillations with the help of a dimension-5 term, 
which leads to a Majorana mass matrix for the neutrinos. By diagonalising this matrix one obtains the 
physical neutrino masses and the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix encoding the 
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strength of the flavour transitions HJIU. The dimension-5 mass term is naturally generated if right- 
handed neutrino fields are added to the Standard-Model (SM) Lagrangian: being gauge singlets these 
fields permit fundamental Majorana mass terms, which are not protected by the SM gauge symmetry 
and can consequently be very large. Integrating out the heavy right-handed neutrinos generates the 
dimension-5 term and the desired small neutrino masses m,^. through the famous see-saw formula 
|I3]^71. With this set-up lepton-flavour violating (LFV) decays of charged leptons like ij — )• ii^ (where 
^1,2,3 = e.,ix,T) occur at unobservably small rates, because the transition amplitudes are suppressed by 
a factor of {ml. — m^J/M^^. This situation is dramatically different in the Minimal Supersymmetric 
Standard Model (MSSM), which contains new sources of flavour violation in the soft supersymmetry- 
breaking sector. To study LFV effects one commonly adopts a weak basis of the (s)lepton multiplets 
in which the lepton Yukawa couplings are flavour-diagonal. The off-diagonal elements of the charged 
slepton mass matrix, Am^y with i,j = 1,2,3 and X,Y = L,R, give rise to LFV decays of charged 
leptons through loops containing a slepton and a neutralino. By confronting MSSM predictions with 
experimental upper bounds on LFV decay rates one can derive constraints on these elements, which 
are usually quoted for the dimensionless parameters 



4y = . (1) 



Here mij^ denotes the i-th diagonal element of the slepton mass matrix with chirality X = L,R. If 
all 5^^Y small, essentially coincides with the corresponding physical i-th generation charged 
slepton mass. A different avenue to constrain the (5^y in Eq. ([T]) are studies of deviations from lepton 
flavour universality (LFU). This approach has been proposed in Ref. El, which exploits the impressive 
experimental precision of 

This notation implies a sum over all three neutrino species. The experimental situation is summarised 
in Tab. [T] The cancellation of the hadronic uncertainties makes the theoretical prediction of Rk very 
clean: Including bremsstrahlung the SM value is given by llOUTTI 



RT = ^P^-^ (1 + '^^QEd) = (2. 477 ± 0. 001) • 10"^ (3) 
mf, (mj^ - mf,y 



The large helicity suppression of the SM contribution to the electronic decay mode makes Rk sensi- 
tive to effects of a charged Higgs boson. In the MSSM the charged-Higgs contribution cancels from 
Rk at tree-level. Yet, as pointed out in Ref. |8|, at the loop level LFU- violating contributions involv- 
ing Am*^y can lead to Rk 7^ R^. It is convenient to parametrise the n — e non-universality in Rk 
in terms of the quantity Ar'*^^ defined as 

Rj^ = i?SM ^ Ar^'-'') . (4) 

Supersymmetric contributions which are linear in J^y cannot interfere with the SM amplitude in 
K — iui, because they lead to a final state with charged lepton and neutrino belonging to different 
fermion generations. Therefore these contributions will necessarily increase r{K — )• iv). In Ref. |i8J a 
mechanism involving the product 5j^^6j^j^ has been proposed to achieve a suppression of r{K eVe) 
and therefore of Rk- Recently, two new observables have been found to constrain the very same 
combination of supersymmetric FCNC parameters fT2l: Firstly, 't Hooft's naturalness criterion has 
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Experiment 


Rk [10-^] 


error ^Rk/rj^ 


PDG 2006 


2. 45 ±0.11 


4.5% 


NA48/2 2003 


2. 416 ±0.043 ±0.024 


2.8% 


NA48/2 2004 


2. 455 ±0.045 ±0.041 


3.5% 


KLOE 


2.55 ±0.05 ±0.05 


3.9% 


Kaon 2007 


2. 457 ±0.032 


1.3% 


PDG 2008 


2. 447 ± 0. 109 


4.5% 


PDG 2010 


2. 493 ± 0. 036 


1.4% 


KLOE 2009 


2. 493 ±0.025 ±0.019 


1.7% 


NA62 Jan 2011 


2. 487 ±0.013 


0.5% 


NA62 Jul 201 1 


2. 488 ±0.010 


0.4% 



Table 1: Experimental values for Rk lfT3l - l2Tll . We use the published result of Ref. JTSl quoted in 
the second-to-last row, which corresponds to Ar'^"*^ = 0. 004 ± 0. 005. The result in the last row, 
reported by NA62 at conferences [20,21 J, corresponds to the full data set collected by NA62 in 2007- 
2008. For the future an experimental accuracy of 6Rk/Rk = 0. 1 — 0. 2% is feasible for the NA62 
experiment ifTTll . 



been applied to the electron mass yielding a non-decoupling upper bound on Secondly, a 

powerful bound on IfJiif^/j/jJ has been derived from the anomalous magnetic moment of the electron. 
The latter constraint decouples, i.e. becomes weaker for larger superpartner masses, but the bounds 
on I'^i'^.'^RijI comparable to the ones found from the electron mass for typical sparticle spectra. In 
Sect. |2] we use the results of Ref. [12] to assess the possible maximal effect on Rk from loop diagrams 
involving Sj^x^^jfj^- In Sect. [3] we study MSSM contributions involving a single power of either 6j^^ or 
S^r>. These contributions feed K evr and therefore increase Rk- In Sect.|4]we conclude. 
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In a two-Higgs-doublet model of type II the following Hamiltonian describes leptonic Kaon decays 
Il23: 



yielding the decay rate 



n7^(l-7')smM(l-7'K 



Ml 



■tan2/3u(l ±7^)51^(1 -75)£ 



T{K ^ lu,) = ^mjmKfk\Vus\'' ( 1 



87r 



m 



K 



2 r 



W^K 



Ml 



(5) 



(6) 



where the second contribution in the square brackets stems from the additional charged-Higgs ex- 
change. At tree level the relative Higgs contribution to the decay rate is independent of the lepton 
flavour and thus cancels in the ratio Rk defined in Eq. However, SUSY loop corrections can 
introduce a dependence on the lepton flavour: In the large tan /3 regime of the MSSM the relation be- 
tween the Yukawa couplings and the measured fermion masses can change significantly, with the loop 
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suppression compensated by a factor of tan (3 ~ 50. In the decoupling limit Msusy ^ v,Mh these 
enhanced corrections arise in a very intuitive way from a loop-induced non-holomorphic Higgs cou- 
pling |[23l . (Here tan /3 = Vu/vd denotes the ratio of the two Higgs vevs, v = \J + v'j = 174GeV 
is the electroweak scale, Mh represents the charged-Higgs boson mass, and Msusy is the mass scale 
of the supersymmetric particles entering the loop diagrams.) In our case of lepton Yukawa cou- 
plings the applicability of the decoupling limit is not clear a priori, because some of the superpart- 
ners involved (e.g. neutralinos) can easily have masses around or even below v. To cover the case 
Msusy ~ v,Mh one must resort to a diagrammatic resummation of tan /3-enhanced corrections, 
which has been worked out for quarks in Refs. II241I25II and for leptons in Refs. 1.12. 261. The desired 
all-order relation between the Yukawa coupling yi and the physical lepton mass rrii is 



(V) 



where is the piece of the one-loop self-energy proportional to m£ tan /3. In the on-shell renormali- 
sation scheme the mass counterterm is just 6m£ = up to terms which are not enhanced by a factor 
of tan /3. Eq. (|7]) is conventionally written as 



. ^ in£ tan f3 

v[l + €e tan p) 



§2 m£ tan /3 
V2Mw 1 + e^tan/3' 



ei tan /3 



me 



(8) 



—yi sin /3 is just the Higgs coupling to right-handed down-type leptons contributing to the second 
terms in the square brackets in Eqs. ^ and Putting everything together one finds 



Rk = rT 



1 — m? 



1 



tan2/3 

^ Mjj (1 + estan/3)(l + eetan/3) 



1 — m? 



1 



K 



Mfj (l + estan/3)(l + e^tan/3)_ 



(9) 



where is the analogue of for the strange Yukawa coupling. This reads 



2m\ tan^ /3 
M'jj\l + es tan/3| 



|l + eetan/3| |l + e^tan/3| 



(10) 



in terms of the notation of Eq. (01). Lepton universality is violated for eg / e^. In the MSSM with 
minimal flavour violation (MFV) the only source of eg 7^ are different values of the selectron and 
smuon masses. A sizable slepton mass splitting between the first and second generation is theoretically 
hard to justify and we do not consider this possibility any further. 

An a priori sizable source of lepton non-universality are the diagrams involving a double insertion 
of LEV mass insertions |8|, see Eig. Instead of the tree level coupling which is proportional to 
rUe tan /? the diagram gives a contribution proportional to nir tan /3. This is the dominant contribu- 
tion to the LEG self energy through double LEV and is given as (with the notation and conventions 
explained in Appendices A, B of Ref. llT2lfl 



Ql 



13 F,(Mf,m? ,m? ,m? ,m? ). (11) 



'The SU(2) partner diagram of Fig.[T]involving a charged Higgs boson was shown in Ref. l8l . 
^Interchanging the chiralities of the two f yields the analogous expression with 5]^''r5'^b.- 
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Figure 1 : Dominant double LFV contribution to the electron mass renormalisation. 



Throughout this paper we choose /x, the gaugino mass parameters, the trilinear SUSY breaking terms, 
and all off-diagonal slepton mass matrix elements real. Therefore e^, are real as well. Eq. (ITTI ) de- 
scribes a non-decoupling effect, because Fq oc M^^y- The contribution of Sf ^ to the resummation 
formulae in Eq. dD is [,12,1 : 

oo fmp + TE^) tan B 

-V smB = — -^^ — ^ - (12) 

V2Mw l + eetan/3 



Writing 



^FV — "T-r tan/3 
1 + er tan /3 



^ = /:"T:. Ai^ (13) 



to adopt a notation similar to Ref. HI the charged-Higgs coupling to the electron changes to 



ig2 rUe 



yesmp = ^j^ tan/3 1 + ^— (14) 

V2Mw 1 + ee tan p \ nie J 

m me « , tan/3 \ 

■tan^ ( 1 + — ^— ^A|^ . (15) 



-v/2Mv(/ 1 + fie tan /3 \ me 1 + e,- tan /3 

Recalling the discussion after Eq. ^ we see that there is a one-to-one correspondence between the 
electron mass counterterm Sf ^ and the enhanced Yukawa coupling in Eq. (fT4l) . In Ref. fV2\ 't Hooft's 
naturalness criterion has been applied to to the electron mass to derive a bound on Within 
theories with MFV the smallness of me is justified by the chiral symmetry gained in the limit ye — )■ 0. 
In our case a second symmetry-breaking parameter, ^r'^ij^'^Jyj is present, and the naturalness principle 
forbids large accidental cancellations between the two contributions to m-e. Demanding < rUe 

one finds the bounds in Tab. [2] 1 12 ], which can be summarised as 

|^iXij|<0.1. (16) 

The authors of Ref. lEl have used A^^ = 0(10^^) to bring Rk into better agreement with the 
NA48/2 result of 2004, see Tab. [T] This value implies a 2000% change in the electron mass which is 
incompatible with the naturalness principle. 

After considering the MSSM we turn our argument into a model-independent analysis: For the nat- 
uralness bound on ye it is inessential in which theory the self-energy Sf ^ is calculated. Any theory 
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with a tree-level Higgs sector corresponding to a type-II 2HDM and a self-energy contribution ^ 
not proportional to affects the charged-Higgs coupling through the finite counterterm 

6y. = ^5m. = "f^"/ ^ ^Sf^. (17) 

rrie rrie 1 + Ce tan p rUe 

If we take the upper bound from the fine tuning argument \6me\ = rUe, the allowed range for y^, lies 
between and twice the tree-level value me/f^- The largest allowed value for — Ar'^"*^ in Eq. (01) 
therefore corresponds to 

Ar'"''' - -A ""^ nS) 



This bound assumes that (as in the MSSM) the muon Yukawa coupling is not substantially affected. 
Taking tan /3 = 50, tan /3 = 0. 3 (which corresponds to a typical loop suppression factor es = 
and tan/3 = 50), and a charged-Higgs mass of Mh = 300 GeV we find Ar^jj^'^p^ = —5 • 10~^, 
which can be probed by NA62 but is not in the 5(t discovery reach of this experiment. Of course 
our consideration equally applies to positive values of Ar^^^LPC i-^- naturalness bound implies 
I Ar^p^l < 5 • 10^^. We have discussed negative contributions, because Ar^"*^ < is an unambigous 
sign of an LFC mechanism, while Ar'^~'^ > can be more easily accomodated with LFV new 
physics as analysed in Sect. [3] In Ref. fTT] (which is an update of Ref. [8|) a thorough analysis 
of several observables in quark and lepton flavour physics has been performed. While most of the 
points in the scatter plots of that paper satisfy the constraint from Eq. ([TSl l. an inclusion of Eq. (fT6l ) 
into the analysis would eliminate the outliers in these plots. Further the use of Eq. ([T6l ) would make the 
results of Ref. [27] less dependent on the anomalous magnetic moment of the muon, whose theoretical 
prediction in the SM involves uncertainties which are not fully understood. 

We note that |5)^^(5jj^^| can also be bounded in a completely different way: The anomalous magnetic 
moment of the electron gives essentially the same bound as Eq. ([T6l ) for MsusY = 500 GeV and 
involves the same supersymmetric particles in the loop as Sf ^ as shown in Ref. lfT2l . Thus relaxing 
the naturalness bound |Sf ^| < rrie to lower Ar^jj^^^pc ^1- CS]> requires the choice of larger 
bino or selectron masses to comply with the electron magnetic moment. Should future NA62 data 
point towards Ar^^*^ < 0, an analysis in conjunction with the electron magnetic moment will place 
correlated lower bounds on these sparticle masses. 

While NA62 is gaining statistics in the forthcoming years, we may expect increasingly better informa- 
tion on Mh and tan /3 from LHC experiments, so that the bound ^r'^^^^^ = — 5 • 10^^ quoted after 
Eq. (fTSl ) may eventually become tighter. Any future NA62 measurement of Ar'^^'^ below Ar^r'j^p^ 
will then establish a more exotic new physics explanation than type-II charged-Higgs exchange, such 
as t-channel leptoquark exchange. 



3 Lepton-flavour violating loop corrections 

Flavour- violating self-energies in the charged-lepton line can induce the decays K — lu^i with i ^ I' 
m. A sizable effect is only possible m K ^ evr, so that LFV self-energies can only increase Rk- 
In this section we estimate the maximal effect of lepton-flavour violating loop corrections to Rk- A 
large correction to Rk involves large tl-tr mixing. We cannot rely on the expansion in ^^/Mgygy 
adopted in Ref. ||8l|27l in this region of the MSSM parameter space, because the f mixing angle Or 
vanishes for f /MsusY — ^ 0. We use the exact formulae of Ref. fTT], which express Ar^^*^ in terms 
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scenario 


.X = 0. 3 


X = 1 


X = 1.5 


X = 3. 


for 


1 


Ml = M2 = mi = mn 


0.261 
0.234 


0.073 
0.059 


0.050 
0.040 


0.026 
0.023 


S\i^„Sr^ > 
^Ir^^^l < 


2 


3Mi = M2 = mi = nijf 


0.301 
0.269 


0.083 
0.067 


0.057 
0.045 


0.029 
0.024 


< 


3 


Ml = M2 = 3m£ = 771 B, 


0.292 
0.235 


0.082 
0.067 


0.057 
0.042 


0.031 
0.027 


^ffl^ii < 


4 


Ml = M2 = If^ = mfl 


0.734 
0.702 


0.210 
0.190 


0.142 
0.127 


0.071 
0.064 


3^RS\:i > 

s]!rS}1 < 


5 


3Mi = M2 = m,i = 3mji 


0.731 
0.693 


0.205 
0.179 


0.137 
0.116 


0.067 
0.054 


> 

< 



Table 2: Different mass scenarios and the corresponding upper bounds for rriR^L denotes 

the average right and left-handed slepton mass, respectively, Mi and M2 the bino and wino masses 
and X = fi/mji. In all scenarios, tan /3 = 50 and sgn /i = +1. 



of the masses of the physical stau eigenstates fi 2 rather than the diagonal elements m?^ ^ of the stau 
mass matrix. Since NA62 runs concurrently with the LHC, it is anyway useful to express Ar^"'^ in 
terms of the physical quantities probed in high-py physics. 

Ti£.^_i.j^ is relevant for the M/^-coupling to leptons if j > i (and thus contributes to thi^eshold correc- 
tions of the PMNS matrix as studied in Ref. [12 ,28 1), whereas for j < i S^^^.^-^ is responsible for the 
correction of the charged-Higgs coupling. The charged-Higgs couplings to leptons including 

an analytic resummation of tan /3-enhanced corrections are listed in Eqs. (31a-c), (32a-c) and (33a- 
c) of Ref. |[T2l . For the decoupling limit Msusy ^ v these charged-Higgs couplings were derived 
earlier in Eqs. (92-95) of Ref. [29], which further uses the iterative procedure of Ref. f30l to resum 
the tan /3-enhanced terms. In this paper we are interested in the tl — )• e/j self-energy contributing to 
Rk as shown in Fig.[2l aiming at constraints on 5^j^^ to be obtained from future NA62 data. With the 
LFV Higgs couplings at hand one can calculate the decay rates summing over all neutrino species and 
compute the ratio. The H~^eVr vertex, which involves the enhancement factor of mr/nie, is the only 
relevant contribution to Rk'- Corrections to the muonic decay mode are irrelevant due to the much 
smaller enhancement factor of mr/m^. In view of the result of the previous section we can further 
rule out large effects in F^^. In the LFV case the deviation from the SM is essentially given as: 



^'LFV 



771 



K 



tan^ 13 



mz 



Mjj{l + es tan/3)^(l + e, 



.tan/3)^ ml 



(19) 



Here S^^_g^ denotes the sum of the two self-energies appearing in Fig.[2l For the tl — )• e/j-transition 
only the two diagrams in Fig. |2]are not suppressed with ye- We explicitly account for stau mixing 
with mixing angle Or € [— 7r/4, 7r/4] and stau mass eigenstates fi = ficosOr — fRsm.6r, f2 = 
TLsmOj- + fjicos Or- We use the following identities (with 777^'^^ = 771t-/(1 + ertan/?) and the 



8 




Figure 2: Dominant LFV self energies in external legs leading to flavour non-diagonal charged-Higgs 
couplings. B and H denote the relevant bino and higgsino components of the neutralinos Xi ■ ■ ■ X4> 
respectively. 



trilinear yl-term set to zero since it is not tan /3 enhanced): 



2Amfj^ _ -2mV Vtan/3 

7 

1 



sin(20O = = "'7 ^"7^ , cos(20,) = (20) 



= 2 V'^^L + ± V (""^i " ""^fi)^ + 4|Am33j|2j (21) 
sgn (m?^ - m|J = sgn (m|^ - m?^) (22) 
The correlation of sgn — m?^ ) with sgn /i and sgn 6*,- can be read off from 

2fim^^ tan /3 
sm(2fc'T-j 

Since we want to use the mixing angle and the mass of the lightest stau, mf,, as inputs we express the 

n2 



masses of the left- and right-handed stau in terms of m? and 6^-: 

2 2 tJ-m^^ tan /3 
sin(26l^) 

= + ^"^ilaXy ^ '^"^^"^ cos(20O) • (25) 



The phenomenologically interesting large values of Ar^^"^ involve large values of Varying 
to larger values with m?^ ^ fixed increases the mass splitting between the two stau mass eigenstates 
fi and f/j and will eventually lower the smaller physical stau mass nif^ below its experimental lower 
bound. We avoid this problem by varying the parameters for fixed mf^ . Treating the flavour violating 
off-diagonal elements up to linear order and including stau mixing we get for the left diagram in 
Fig. 12 

^r^-e^ = - ^fn ^RR Or COS Qr (/i (M^, m?^, m?J - /i ( Af f , rn?^ , ) 

«! ,^ t-is mT-tan/3 „ /,,2 2 2 2 \ 

= _ _ Ml mg^ m^^ /X ^^^^^^^^ h [Ml , m^^ , m^^ , m^J . 

Here and in the following we need the loop functions 

^ xylnf+xzlnf +j/zln| fiix,y,z) - hix,y,w) 

fi{x,y,z) = — , f2{x,y,z,w) = . (27) 

[X — y){x — z){y — zj z — w 
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For the right diagram in Fig. [2] we get the following contribution: 

TL en 1 + er tan ^ ^ ^ (28) 

• [sin2 0^/2 {Ml i?,m\ , + cos^ ^./a (Mf, m?^ , m?^)] . 

In order not to get negative slepton masses, must for sure be smaller than 1. Here, we only 

consider a single flavour-violating mass insertions. Contributions with double mass insertion, e.g. 
^f^L^^RR ^^"^ relevant for /i — )■ 67 IBT]. In principle, T?^^_^^ = + is sensitive to the 

RR-element. However, the relative minus sign is the origin of a possible cancellation in certain region 
of the parameter space. In this approximation the sensitivity to 5^^^ vanishes if jj? = m?^ cos^ Or + 
sin^ Or- In the case with = the cancellation occurs for ij? = mf^. This feature was already 
discovered in Ref. {3T]. In Refs. fS^.Tf] the decoupling limit AfsusY ^ v is adopted. In this limit 
fL,R appear in the loop functions instead of fi^2 and the tl-tr flip in the bino diagram is incorporated 
within the mass insertion approximation (MIA). This results in a simplified version of our Eqs. (l26l ) 
and (l28]l; in Eq. the square bracket simphfies to /2 [Mj, ,m1^,ml^ Since Or vanishes for 
-^SUSY S> V, the consideration of large stau mixing requires to go beyond the decoupling limit and 
beyond MIA. Furthermore the stau and bino masses can still be smaller than v; in fact the interesting 
region of the parameter space probed by NA62 comes with light bino and staus. 

We now estimate the maximal allowed LFV effect in Rk including stau mixing, which depends very 
much on /itan/3. At the end of Sect. |2] we already concluded that effective LFC effects are typically 
below the experimental sensitivity. In our plots and numerical examples we use the following values 
for the smaller stau mass ■mf^ and the bino mass parameter Mi. 

rrif, = 120 GeV, Mi = 100 GeV (29) 

These values are consistent with the experimental lower bounds of 46 GeV for the neutralino masses 
and 81. 9 GeV for nif^ |[T3l . The heavier stau mass m^^ is then calculated from the mixing angle Or 
and /i. For the off-diagonal element Am^^j^ = me^mf^5)j^^ we choose for simplification mg^ = 
200 GeV, with m-f^ also calculated from Or, /x, and mf^. With this choice the bino diagram increases 
with /X, since m^^ and mf^ increase too. (Setting instead AmJ^j = would lead to a finite 

limit for S^_g^ for fi — )• 00.) Furthermore, we choose tan/3 = 50 and /j, to be real and positive. 
With this chosen input parameters we can analyse the dependence of Ar^~^ on Or, /U and Si^rr- For 
small values of /i the higgsino-bino diagram dominates, but the maximal value is rather small. For 
large /i (and other SUSY masses fixed) the higgsino-bino diagram tends to zero whereas the pure 
bino-diagram can become sizeable. Without stau mixing this diagram would not contribute at all. 
Thus, in order to get any sizeable effect, especially for large values of /i, one has to take stau mixing 
into account. With this setup the largest effect in T,^^_^^/mr comes with a relatively large mixing 
angle of Or ~ 26°. In case of the higgsino-bino diagram, stau mixing is not important, this diagram is 
maximal for small /x, but nevertheless approximately one order of magnitude smaller than the largest 
values found for the pure bino diagram.For Or ~ 26° and /x > one has fi = fi, which moreover is 
dominantly left-handed. 

The maximal possible deviation of Rk from the SM prediction is visualized in Fig. [3] where Ar'^~'^ 
is plotted for 6j^j^ = 0. 5, Mh = 500 GeV and tan /5 = 50 as a function of Or and /x using typical 



'The relation between the notation in Ref. I8I27I and ours is Et^, -er = j^f^j^A^jf. E.g. the second term in Aff cor- 
responds to the bino diagram with the following simplification compared to Eq. i26i : sin 6t cos dr (/i {M^ , rrig^ , rrii^ ) — 
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Figure 3: Ar^-^ for 5^^j^ = 0. 5, Mh = 500 GeV and tan/3 = 50. Left: As a function of 9r for 
different values of /i: 800 GeV (red), 400 GeV (blue dashed), 200 GeV (green dotted). Right: In 
dependece of /i for different values of Or'. 26° (red), 45° (blue dashed), —18° (green dotted). 



values of tan ^5 = 0. 3 and e^, tan /3 = e,- tan (3 = —0. 07. The discontinuity at 0^ = just comes 
from the fact that fi and f2 change their roles as heavier and lighter staus. In order to find Ar'^"'^ for 
different values of (5)yj, Mh and tan /3 one must rescale those plots using that Ar'^~'^ is quadratic in 
6]ljl, oc and oc tan^ /3. One gets a maximal effect of 0.6% for our chosen point of fi = 800 GeV, 
Mh = 500 GeV, (5jf^ = 0. 5, 9^ = 26° and tan /3 = 50, which is already in the reach of NA62. In 
the range of 500 GeV < fi < 900 GeV a handy approximate formula (with an error of 7%) for the 
maximal effect (occuring at 6^- = 26°) can be found 



^W,LFV-U-UU6(^ M^, M 50 {o.bj UoOGevJ 



valid for nif, = 120 GeV, Mi = 100 GeV, mg^ = 200 GeV. (30) 



2 



If one varies the lightest stau mass, Ar^ scales approximately as (120 GeV/mfJ in the range 
100 GeV < nifj < 250 GeV. The dependence on Mi is roughly hnear for 50 GeV < Mi < 100 GeV 
and the prefactor in Eq. (l30l ) decreases from 0.006 to 0.0028 if Mi is lowered to 50 GeV. Above 
100 GeV the Mi dependence flattens off with a maximum at 200 GeV, at which the prefactor of our 
approximate formula becomes 0.0078. A full investigation of the dependences of Ar'^"'^ on Mi and 
mg^ requires the use of the exact expression, obtained by adding the quantities in Eqs. ( [28] ) and ( [261 ) 

to find and inserting the result into Ar^~*^ in Eq. ([T9l ). 

In Fig. |4]we show the dependence of Ar^^^ on fi, Mh, tan /3 and 6j^j^ for 6^- = 26°. It is possible 
to reach an effect of 0(0. 5%), whereas for vanishing mixing or small jj, it is hardly possible to reach 
the experimental sensitivity. To derive contraints on 6j^j^ Fig. [5] might be useful. We show regions in 
the (^MH,fJ-, tan/3,(5|yj) parameter space where Ar'^^^ reaches the future experimental sensitivity 
of 0.2%. In Ref. [12] it is pointed out that even large values of tan /3 = 100 are compatible with the 
requirement of pertubative bottom Yukawa coupling. 

In order to estimate the contribution of tl — )• transition to Rk in Ref. |l8|, the authors used 
A|^ = 5 • 10"'' which is related to our notation by Sri-e^ = i+I^t^if/^ ^fi "^tI-cr l^r = 
0. 025. However, such large values correspond to quite specific points in the MSSM parameter space, 
especially extremely large ^. Thus, in order not to overestimate the effect and to avoid trouble with 
too small slepton masses one should rather take at most 



yX" 



/nir « 0. 01 (meaning A^^^^^ 
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A ^''"" = 0.002 A /■'' "=0.002 




200 300 400 500 600 200 300 400 500 600 700 800 900 



Mh [GeVl [GeVl 

Figure 5: For different values of 6]^j^ = 0. 15 (yellow), 0. 25 (red), 0. 5 (green), 0. 75 (blue) (from top 
to bottom) we plot the regions in which Ar^"*^ is below the future experimental sensitivity of 0.002 
in the MH-tan (3 plane with = 800 GeV (left) and in the ^-tan f3 plane with Mh = 500 GeV 
(right) and stau mixing angle O-j- = 26°. I.e. if 6]fj^ = 0. 25, the white and yellow areas correspond to 
Art"-" > 0. 002. 
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2 • 10^^). Taking this maximal value and further tan/3 = 50 and a charged-Higgs mass of Mh = 
500 GeV we end up with Ar^~^'^LFV ~ 0. 007, which is within the experimental sensitivity of the 
NA62 experiment. 

In Ref. lUPTl the 95%CL exclusion region for Rk for three different values of A|j (1 • 10"'^, 5 • 
10"^, 1 • 10"^) is shown in the (A%± , tan /3) region and compared with the constraints from B — )• tv, 
B Xs^, Rf_t23 = ^{K fiv)/T{K TT^fiu) (331 and direct searches. According to the 
discussion in the preceding paragraph we prefer to take A|j^ = 2 • 10~^ as the maximal value, so 
that our excluded (Mj:^±,tan/3) region is smaller than the one in Ref. fT9l. In the following we 
use |l34H37l 

B{B ruf^ = 1. 13 • 10-^ (i^)' (200^)' , (31) 
B{B Tuf''^ = (1. 64 ± 0. 34) • 10"^ , (32) 

B[B ^ rr/)™^^ = 

with eo ~ eg « ertan/? = —0. 07and ji = 800 GeV as above. In Fig. |6]we plot the region 
in the in the Mj:^±-tan/3 plane satisfying /\r^~^ < 0. 5% for three different values of 6j^j^ (0.15, 
0.25, 0.5). Overlaid are the constraints from B — )• tu, R^23 = ^{K — fiv)/T{K — ir^fiu) and 
direct searches. The prediction of B{B — )• tu) within the SM and the MSSM requires the 
knowledge of \ Vub\- Determinations of \ Vub\ from different quantities result in substantially different 
numerical predictions. For discussions of this "14fc-puzzle" see Refs. Il34[|371l38l . For our analysis we 
consider two extreme scenarios: First, in the left plot of Fig. |6] is determined such that the SM 
prediction of B{B — )• tu) is equal to the experimental value. Using fs = (191 it 13) MeV one gets 
\Vub\ = (5. 04 ± 0. 64) • 10-3 [34|. In the plot we set B{B tu)^^ = B{B ri^y'P and use the 
experimental 3cr-region. Second, in the right plot of Fig.|6]l4{, is fixed to the best-fit value of a global 
fit to the unitarity triangle li38l . An essential assumption of the second scenario is the absence of new 
physics in the CP asymmetry Aqp{B — )• J/ipKs), from which the angle (3 of the unitarity triangle is 
determined: Then |Kb| oc |Kblir| leads to = (3. 41 ± 0. 15) • 10"''^ |34|. In this case the SM 
central value B{B — tv)^^ = 0. 75 • 10"^ is much lower than the experimental value and the SUSY 
contribution makes it even smaller. Therefore now B{B — )• tu) is much more constraining than in the 
scenario of the left plot. Using the 2a region of the experimental value one would exclude the whole 
{Mfj±,tan region (except for a very narrow strip with tan /3 0.3Mjj+y, the constraint from 
the 3(T region is shown in the right plot of Fig.[6l In the future also B{B — )• Dtv) |[37l[39 .40l and 
B{B — )• TTTi/) BTl will probe charged-Higgs effects and will eventually shed light on the situation. 

We finally mention two related studies: A prospective error of 0.12% of the NA62 experiment at 
CERN is used in Ref. 1421 . The pai^ameter scan in this paper respects — 0- O-*-' agreement 

with our result. In Ref. f43] it is pointed out, using a general effective theory approach, that in 
models with Minimal Lepton Flavour Violation (MLFV and MFV-GUT) the effects are too small to 
be observed. 

4 Conclusions 

The NA62 experiment has the potential to discover new sources of lepton flavour violation by testing 
lepton flavour universality through a precision measurement of Rk = T {K — )• ev) /T {K — )• ixu) |i8J. 



rriB 



tan^ (3 



mu+ J (1 + eo tan /3)(1 + er tan/3) 



B{B 



TV 



,SM 



(33) 



4 Conclusions 



13 




200 400 600 800 1000 200 400 600 800 1000 

Mh* [GeV] [GeV] 

Figure 6: Regions with Ar^"^ > 0. 5% for 6]^j^ = 0.15 (darkblue), 0.25 (blue and darkblue), 0.5 
(lightblue, blue, and darkblue). Overlaid in red: exluded by LEP ff+ searches; yellow dashed: 
3a exclusion limit from R^23', green: 3cj exclusion limit from B ^ tv (left: using \Vub\ = 
(5. 04 ± 0. 64) • 10-3; right: using \ Vuh\ = (3. 41 ± 0. 15) • lO^^). 



This kind of new physics dominantly affects the decay rate T{K — )• ev). A lepton-flavour conserving 
(LFC) mechanism changing T{K — )• efg) may suppress or enhance Rk, while new lepton-flavour 
violating (LFV) decay modes such as T{K — )• evr) can only enhance Rk over its SM value. In this 
paper we have studied Ar^^^"^ = Rk/R^ - 1 in the MSSM, extending the analyses of Refs. Il8ll27l. 

The LFC contribution to Ar^^"^ is driven by the parameter combination In Ref. | [T2l it has 

been found that upper bounds on | can be derived from naturalness considerations of the 

electron mass and from the precise measurement of the anomalous magnetic moment of the electron. 
(Coincidentally, these two quantities give very similar constraints.) In Sect. [2] we have found that 
these bounds imply | Ar^pf?! < 0. 005 and thereby challenge the large values for | Ar^p^j considered 
in Ref. [8]. At the same time our result is fully compatible with the range for Ar^p^? advocated in 
Ref. ll27l . The naturalness bound extends beyond the MSSM to a larger class of models, namely those 
with the tree-level Higgs sector of a 2HDM of type II. 

The LFV contribution to Ar^""^ can be larger, because a non-zero parameter suffices to open 
the decay channel K — eur and 6jfj^ is only poorly constrained from other processes. We have 
calculated Ar^p^ in Sect. [3] and found that the proper inclusion of f^^-f/j mixing is essential. The 
analytical expressions in Refs. ||8l|23| include the tl — tr flip using the mass insertion approximation 
instead of the exact diagonalisation of the stau mass matrix. The interesting region of parameter space 
probed by NA62 corresponds to large values of ^ and a sizable stau mixing angle Or and in this region 
the left (bino) diagram in Fig.|2]is dominant. The formulae derived by us are also valid beyond the 
decoupling limit MsusY oo, in which 9r vanishes. In order to facilitate the combination of future 
NA62 results with limits or measurements from high-p^ experiments, we have expressed Ar^py in 
terms of the mass mf^ of the lightest stau eigenstate and the mixing angle 6^. For example, for 
tan/3 = 50, /i = 800 GeV, (5|fj^ = 0. 5, a charged-Higgs mass of Mh = 500 GeV, ruf^ = 120 GeV, 
a bino mass of Mi = 100 GeV and a right-handed selectron mass of mg^ = 200 GeV we find a 
maximal value of Ar^py = 0. 006 corresponding to Oj- = 26°. In Eq. (l30l ) we have derived an easy- 
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to-use formula expressing Ar^jS^ in terms of the relevant MSSM parameters. Finally we have plotted 
the regions of the MSSM parameter space probed by Rk and briefly compared the result with the 
constraint from other observables such as B{B — )• tu). 

Acknowledgements 

The presented work is supported by BMBF Grant No.05H09VKF. J.G. acknowledges the financial 
support by Studienstiftung des deutschen Volkes and the DFG cluster of excellence "Origin and Struc- 
ture of the Universe". U.N. is grateful for the hospitality of the Institute for Advanced Study of 
Technische Universitat Miinchen, where this paper has been completed. 

References 

[1] B. Pontecorvo, Inverse beta processes and nonconservation oflepton charge, Sov. Phys. JETP 7 
(1958) 172-173. 

[2] Z. Maki, M. Nakagawa, and S. Sakata, Remarks on the unified model of elementary particles, 
Prog. Theor Phys. 28 (1962) 870-880. 

[3] P. Minkowski, ^ ^ ej at a Rate of One Out of 1 -Billion Muon Decays?, Phys. Lett. B67 (1977) 
421. 

[4] T. Yanagida, Horizontal symmetry and masses of neutrinos, ConfProc. C7902131 (1979) 95. 

[5] S. Glashow, The future of elementary particle physics, NATO Adv. Study Inst.Ser.B Phys. 59 
(1980) 687. Preliminary version given at Colloquium in Honor of A. Visconti, 
Marseille-Luminy Univ., Jul 1979. 

[6] M. Gell-Mann, P. Ramond, and R. Slansky, Complex spinors and unified theories, ConfProc. 
C790927 (1979) 315-321. To be pubhshed in Supergravity, P van Nieuwenhuizen and D.Z. 
Freedman (eds.). North Holland Publ. Co., 1979. 

[7] R. N. Mohapatra and G. Senjanovic, Neutrino Mass and Spontaneous Parity Violation, 
Phys.Rev.Lett. 44 (1980) 912. 

[8] A. Masiero, P. Paradisi, and R. Petronzio, Probing new physics through n — e universality in 
K lu, Phys. Rev D74 (2006) 011701, |hep-ph70 5TT2 8 9| . 

[9] M. Finkemeier, Radiative corrections to 7r;2 and K12 decays, Phys. Lett. B387 (1996) 391-394, 
ihep-ph/9505434| . 

[10] V. Cirigliano and I. Rosell, Two-loop effective theory analysis of pi K — )• ei>e(7) branching 
ratios, Phys.Rev.Lett. 99 (2007) 231801, MarXiv: 0707 . 34391 . 

[1 1] V. Cirigliano and I. Rosell, ir/K ev branching ratios to 0{e'^p'^) in Chiral Perturbation 
Theory, JHEP 10 (2007) 005, Il0707 . 446411 . 



[12] J. Girrbach, S. Mertens, U. Nierste, and S. Wiesenfeldt, Lepton fiavour violation in the MSSM, 
JHEP 05 (2010) 026, llarXiv : 910 . 2 6631 . 



REFERENCES 



15 



[13] Particle Data Group Collaboration, K. Nakamura et. al, Review of particle physics, J. Phys. 
G37 (2010) 075021. 

[14] G. Isidori, KAON 2007: Conference Summary, PoS KAON (2006) 064, 
[larXiv: 70 9 ■24381 . 

[15] KLOE Collaboration, F. Ambrosino et. al. Precise measurement of 

B{K eu{'j))/B{K ^1^(7)) and study of K ei/7, Eur.Phys.J. €64 (2009) 627-636, 
llarXiv: 907 ■ 35941 . 

[16] NA62 Collaboration, E. Goudzovski, Lepton flavour universality test at the CERN NA62 
experiment, Nucl.Phys.Proc.Suppl. 210-211 (2011) 163-168, llarXiv: 10 08 .121911 . 

[17] E. Goudzovski, Talk at NA62 Physics Handbook Workshop, CERN (2009). 

[18] NA62 Collaboration, C. Lazzeroni et. al.. Test of Lepton Flavour Universality in — )■ l^v 
Decays, Phys.Lett. B698 (2011) 105-114, llarXiv: 1101 . 4805L 

[19] E. Goudzovski, Lepton flavour and number violation with K decays at CERN, Talk at 46th 
Rencentres de Moriand, La Thuile, Italy (2011). 

[20] E. Goudzovski, Kaon experiments at CERN: NA48 and A?A62,larXiv : 1112 . 53 65] 

[21] NA48/2 and NA62 Collaboration, E. Goudzovski, Kaon programme at CERN: recent results, 
TrXiv: 1111 .2818. 

[22] W.-S. Hou, Enhanced charged Higgs boson effects in — )■ tu, and h ^ tu + X, Phys. 
Rev D48 (1993) 2342-2344. 

[23] L. J. Hall, R. Rattazzi, and U. Sarid, The Top quark mass in supersymmetric SO(IO) unification, 
Phys. Rev D50 (1994) 7048-7065, [hep-ph /9306309 ]. 

[24] M. S. Carena, D. Garcia, U. Nierste, and C. E. M. Wagner, Effective Lagrangian for thetbH^ 
interaction in the MSSM and charged Higgs phenomenology, Nucl. Phys. B577 (2000) 88-120, 
||hep-ph/ 9 912 5161. 

[25] L. Hofer, U. Nierste, and D. Soberer, Resummation of tan ^-enhanced supersymmetric loop 
corrections beyond the decoupling limit, JHEP 0910 (2009) 081, llarXiv : 907 . 5408 1. 

[26] S. Marchetti, S. Mertens, U. Nierste, and D. Stockinger, Tan (^-enhanced supersymmetric 
corrections to the anomalous magnetic moment of the muon, Phys. Rev. D79 (2009) 013010, 
llarXiv: 80 8 . 15Tol . 

[27] A. Masiero, R Paradisi, and R. Petronzio, Anatomy and Phenomenology of the Lepton Flavor 
Universality in SUSY Theories, JHEP 11 (2008) 042, llarXiv: 0807 . 47211 . 

[28] A. Crivellin and J. Girrbach, Constraining the MSSM sfermion mass matrices with light fermion 
masses, Phys. Rev D81 (2010) 076001, llarXiv: 1002 . 02271 . 

[29] J. Hisano, M. Nagai, and P. Paradisi, Flavor effects on the electric dipole moments in 

supersymmetric theories: A beyond leading order analysis, Phys. Rev. D80 (2009) 095014, 
liarXiv: 0812 . 4283il . 



16 



[30] A. J. Buras, P. H. Chankowski, J. Rosiek, and L. Slawianowska, AM^^s, ^ /^^/" ^""^ 

B ^ Xg^ in supersymmetry at large tan /?, Nucl. Phys. B659 (2003) 3, [ hep-ph/0210145| . 

[31] J. Hisano, M. Nagai, P. Paradisi, and Y. Shimizu, Waiting for // — )• ejfrom the MEG 
experiment, JHEP 0912 (2009) 030, IfarXlv: 0'90 4 . 20 801 . 

[32] P. Paradisi, Higgs-mediated e ^ p. transitions in II Higgs doublet model and supersymmetry, 
JHEP 08 (2006) 047, |hep-ph/0601100i . 

[33] M. Antonelli, V. Cirigliano, G. Isidori, F. Mescia, M. Moulson, et. al.,An Evaluation of\Vus\ 
and precise tests of the Standard Model from world data on leptonic and semileptonic kaon 
decays, Eur.Phys.J. C69 (2010) 399-424, [ arXiv: 1005 . 2323J . 

[34] U. Nierste, Flavour Physics, supersymmetry and GUTs, talk at The Role of Heavy Fermions in 
Fundamental Physics, April 1 1 -14, Portoroz, Slovenia (2011). 

[35] U. Nierste, Flavour physics, supersymmetry and grand unification, talk at 46th Rencontres de 
Moriond on Electroweak Interactions and unified Theories, 13-20 Mar 2011, La Thuile, Italy 
(2011) larXi v: 1107 . 06211. 

[36] BABAR Collaboration, R. Barlow, Experimental status of B — tv and B — £i'{'y), 
arXiy ;11.02 j_12.67. 

[37] U. Nierste, S. Trine, and S. Westhoff, Charged-Higgs effects in a new B — Dtu differential 
decay distribution, Phys.Rev D78 (2008) 015006, [arXiv: 08 01 . 4 938]. 

[38] A. Lenz, U. Nierste, J. Charles, S. Descotes-Genon, A. Jantsch, et. al, Anatomy of New Physics 
inB - B mixing, Phys.Rev D83 (2011) 036004, |arX i,v : IJjO 8 , 15 93]. 

[39] J. F. Kamenik and F. Mescia, B — )• Dtu Branching Ratios: Opportunity for Lattice QCD and 
Hadron Colliders, Phys.Rev D78 (2008) 014003, [arXiv: 0802. 3790]. 

[40] M. Tanaka and R. Watanabe, Tau longitudinal polarization in B — )• Dtu and its role in the 
search for charged Higgs boson, Phys.Rev. D82 (2010) 034027, l arXiv: 10 05 .4 306]. 

[41] A. Khodjamirian, T. Mannel, N. Offen, and Y.-M. Wang, B Tiivi Width and \ Vub\from QCD 
Light-Cone Sum Rules, Phys.Rev D83 (2011) 094031, l arXiv: 110 3 . 2 6551- 

[42] J. Ellis, S. Lola, and M. Raidal, Supersymmetric Grand Unification and Lepton Universality in 
K lu Decays, Nucl.Phys. B812 (2009) 128-143, r arXiv: 80 9 . 521 1|. 

[43] A. Filipuzzi and G. Isidori, Violations of lepton-fiavour universality in P ^ Iv decays: a 
model-independent analysis, Eur Phys. J. C64 (2009) 55-62, la rXiv: 0906.30241 . 



